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The antiferromagnetic Heisenberg model on the dodecahedron possesses a number of ground state 
magnetization discontinuities in a field at the classical and quantum level, even though it lacks mag¬ 
netic anisotropy. Here the model is considered for two dodecahedra coupled antiferromagnetically 
along one of their faces, as a first step to determine the magnetic response of collections of fullerene 
molecules. The magnetic response Is determined from the competition among the intra-, interdo- 
decahedral exchange and magnetic field energies. At the classical level the discontinuities of the 
isolated dodecahedron are renormalized by the interdodecahedral coupling, while new ones show 
up, with the maximum number of ground state discontinuities being six for a specific range of the 
coupling. In the full quantum limit where the individual spin magnitude s = §, there are two 
ground state discontinuities originating in the single discontinuity of the isolated dodecahedron, and 
another one due to the Intermolecular coupling, generating a total of three discontinuities which 
come one right after the other. These results show that the magnetic response of more than one 
dodecahedra interacting together is quite richer than the one of a single dodecahedron. 

PACS numbers: 75.50.Xx Molecular magnets, 75.10.Jm Quantized spin models, including quantum spin 
frustration, 75.10.Hk Classical spin models, 75.50.Ee Antiferromagnetics. 


I. INTRODUCTION 

The antiferromagnetic (AFM) Heisenberg model 
(AHM) has been extensively investigated in the recent 
decades as a prototype for strongly correlated electronic 
behavior [I|, Special attention has been reserved 

for lattices and clusters with frustrated connectivity, 
which in combination with low spatial dimensionality 
and strong quantum fluctuations can lead to unexpected 
magnetic behavior This includes phases without 

conventional order, such as the spin-liquid phase, non¬ 
magnetic excitations inside the singlet-triplet gap, and 
magnetization plateaux and discontinuities in the re¬ 
sponse to an external magnetic field. 

In the case of lattices frustration manifests itself in the 
form of magnetization plateaus and discontinuities 

Finite clusters can also exhibit magnetization dis¬ 
continuities. A class of molecules associated with mag¬ 
netic frustration when the AHM is considered for spins 
sitting on their vertices are the fullerene molecules. These 
are hollow carbon molecules that come in the form of 
closed cages 0, with structures that can possess high 
spatial symmetry. They are made of 12 pentagons and 
a number of hexagons which varies with the numbers of 
vertices n as ^ — 10. Frustration originates in the pen¬ 
tagons and decreases on the average with n. The poly¬ 
gons that make up the molecules share their edges, while 
each vertex is three-fold coordinated. Maybe the most 
representative member of the class is Cgo, which has the 
shape of a truncated icosahedron and the spatial symme¬ 
try of the largest point symmetry group, the icosahedral 
group Ifi- Ceo was found to superconduct when doped 
with alkali metals [l^, and lies in the intermediate U 
regime of the Hubbard model [ 30 . It was shown that 
in the large U limit of the Hubbard model, the AHM, 
the particular connectivity of Ceo leads to a disconti¬ 


nuity of the magnetization as a function of an external 
magnetic field in the classical ground state 0. This is 
particularly appealing, as the model lacks any magnetic 
anisotropy. A classical ground state magnetization dis¬ 
continuity was also found for the dodecahedron, which 
is the smallest member of the fullerene class and has 20 
vertices and also Ih spatial symmetry. The investiga¬ 
tion was then extended more generally to molecules of 
Ih symmetry. First, it was shown that the dodecahedron 
has in fact a total of three ground state magnetization 
discontinuities in a field at the classical level, and one and 
two ground state discontinuities respectively at the full 
quantum limit of individual spin magnitude s = i and 1 
[3 [3| • It was also shown that the total number of clas¬ 
sical ground state magnetization discontinuities for the 
truncated icosahedron is in fact not one but two, and this 
is a general feature of fullerene molecules of Ih symmetry. 
It was established that another general feature of the Ih 
fullerenes is the high-field discontinuity for s = For 
relatively small fullerene clusters of different symmetry 
only pronounced magnetization plateaus were found for 
s = 1/2 [0. In addition the icosahedron, which is not a 
member of the fullerene family but is the smallest cluster 
with Ih symmetry, has a classical magnetization discon¬ 
tinuity in its lowest energy configuration which persists 
for lower values of s 0, [3] . It must also be noted that 
the on-site repulsion has been found to be stronger for 
C 20 than Ceo in numerical calculations [U, [01 1 provid¬ 
ing further support for the validity of the AHM as a very 
good approximation of the Hubbard model for the do¬ 
decahedron. Work closely related to the above has also 
been published in the literature [3 - [0 . 

While Ceo spontaneously forms in condensation or 
cluster annealing processes [0, this was not the case 
for C 20 , which was eventually produced in the gas phase 
by Prinzbach et al. [0. C 20 has been synthesized in 
the solid phase by Wang et al., who produced a hexago- 
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nal closed-packed crystal [^, while Iqbal et al. synthe¬ 
sized an fee lattice with C 20 molecules interconnected by 
two additional carbon atoms per unit cell in the intersti¬ 
tial tetrahedral sites [s^- On the other hand quasi one¬ 
dimensional structures have been realized experimentally 
only a few times for Ceo molecules, due to their highly 
anisotropic configuration. This included peapods where 
Ceo molecules were introduced in carbon nanotubes, one¬ 
dimensional Ceo structures aligned along step edges on 
vicinal surfaces of metal single crystals, and chains of Ceo 
on self-assembled molecular layers [331 - 1^ . Most recently 
Ceo molecules were arranged on chain structures of width 
two to three molecules on rippled graphene [s^- In ad¬ 
dition, as few as two Ceo molecules have been considered 
to link and form a dumbbell structure [s^. 

The formation of J/j-fullerene lattice structures poses 
the question of the influence of intermolecular interac¬ 
tions on isolated molecule properties. Considering inter¬ 
actions again at the level of the AHM, it is of interest 
to determine if the appealing ground state magnetic fea¬ 
tures of a single dodecahedron survive in a lattice-type 
setting, and if the addition of intermolecular magnetic 
exchange introduces extra features in the magnetization. 
This is the main question undertaken in this paper. More 
specifically the case of two dodecahedra is investigated, 
with intramolecular interactions exactly as in the isolated 
dodecahedron case, while the two dodecahedra are con¬ 
nected along with one of their faces with a varying ex¬ 
change interaction. The properties of the ground state of 
the whole cluster are mapped as a function of the inter¬ 
molecular exchange constant and an external magnetic 
field. For weak intermolecular coupling the response is 
mainly determined by the isolated dodecahedra, while 
for strong by the dimer-type interaction between spins 
belonging to different dodecahedra. 

For relatively small intermolecular coupling the three 
classical discontinuities of the lowest energy configuration 
of the isolated dodecahedron survive. Simultaneously a 
new low-field magnetization discontinuity appears, which 
relates to the AFM coupling of the two molecules and its 
competition with the magnetic field. This discontinuity 
survives up to the dimer limit. For stronger coupling even 
more discontinuities appear, producing a rich structure 
for the classical magnetic response of the two dodeca¬ 
hedra system. For a specific range of interdodecahedral 
coupling the total number of ground state discontinuities 
goes up to six. The spins associated with the interdo¬ 
decahedral interaction do not necessarily increase their 
projection along the field as the latter increases, due to 
their unfrustrated AFM interaction. Finally, one of the 
discontinuities becomes one of the susceptibility close to 
the dimer limit. 

In the s = ^ case the isolated dodecahedron disconti¬ 
nuity generates two ground state discontinuities for the 
coupled dodecahedra. In addition, a third discontinu¬ 
ity appears due to the interdodecahedral coupling. All 
three discontinuities persist for smaller values of the in¬ 
terdodecahedral coupling and appear one right after the 


other. This generates a ground state magnetic response 
which for a considerable range of magnetic fields is as¬ 
sociated with magnetization discontinuities of the total 
spin along the z axis equal to AS^ = 2, instead of the 
typical = 1 for quantum mechanics. 

The plan of this paper is as follows: In Sec. |n]the AHM 
for the system of the two dodecahedra is introduced, and 
in Sec. imi its lowest energy configuration is calculated 
for classical spins. Section m considers the s = ^ case, 
with perturbation theory for weak and Lanezos diago- 
nalization for arbitrary values of the interdodecahedral 
coupling. Finally Sec. IVl presents the conclusions. 


II. MODEL 

The AHM for two linked dodecahedra (Fig. [IJ is: 

H = Si ■ Sj + S20-|-i • ^20-l-j) 

<ij> <20+i,20+j> 

5 N 

+ j'^Si ■ S20+i - h^Si (1) 

i=l i=l 

The total number of spins for the two dodecahedra is 
TV = 40, with the first dodecahedron containing spins I 
to 20 and the second 21 to 40. The first two sums run 
over nearest-neighbor spins within the same dodecahe¬ 
dron, with i and j running from 1 to 20. The second 
sum connects the spins on the faces of the two dodeca¬ 
hedra that are taken to be directly opposite each other, 
with i an index counting the spins on these two faces. 
The magnetic field h is taken to be directed along the z 
axis. The system interpolates between two independent 
dodecahedra for J' — 0, and five independent dimers for 
J = 0. The ratio a = -j^rjr is defined, which correspond¬ 
ingly varies between 0 and 1. 


III. CLASSICAL SPINS 

First the spins of Hamiltonian (ED are taken to be clas¬ 
sical d^. When J' — Q the ground state magnetization 
of each isolated dodecahedron has three discontinuities 
in the field, occuring when = 0.26350, 0.26983, and 

^sat 

0.73428, with = (3 -I- \/5) J the saturation field of an 
isolated dodecahedron 0. The symmetry of the ground 
state configuration does not necessarily increase with 
the magnetic field. In zero field nearest-neighbor spins 
are not antiparallel due to frustration, and the nearest- 
neighbor correlation in each dodecahedron equals — 

On the other hand, if the J' bonds between different do¬ 
decahedra were to be considered alone their spins would 
be antiparallel in the lowest energy state, consequently 
the non-frustrated dimer bonds should be less suscepti¬ 
ble to an external field in comparison with the frustrated 
intradodecahedral bonds. 
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For finite J' and zero field the relative spin orienta¬ 
tions in each dodecahedron in the ground state do not 
change with respect to the noninteracting case, while 
spins connected via J' bonds align themselves in an an¬ 
tiparallel fashion and the energy is —20\/5J— 5J'. Once 
the magnetic field is switched on the competition among 
the intra-, interdodecahedral exchange and magnetic field 
energies determines the lowest energy configuration. The 
magnetization discontinuities associated with an isolated 
dodecahedron survive the interdodecahedral coupling, 
while new ones emerge. The location of all the lowest en¬ 
ergy configuration discontinuities with respect to a and 
is shown in Fig. [2] (hsat is the saturation field of 

the two dodecahedra, which is a function of ^). Apart 
from a close to 1, where the second discontinuity with re¬ 
spect to the field strength becomes the sole susceptibility 
discontinuity, there are never less than four magnetiza¬ 
tion discontinuities, showing that the introduction of the 
interaction between the dodecahedra enriches the mag¬ 
netic response for any coupling strength. The maximum 
number of magnetization discontinuities occurs for a ^ | 
and is equal to six. The discontinuities occur mostly for 
decreasing with increasing a, until they eventually 
disappear in the dimer limit a = 1. The inaccessible 
magnetizations per spin which fall between the edges of 
each discontinuity are plotted in Fig. [H The correspond¬ 
ing magnitudes of the magnetization change per spin are 
shown in Fig. 0] The width of the nonaccessible magne¬ 
tizations is not necessarily monotonic with a. 

Once J' becomes non-zero, apart from the three dis¬ 
continuities of J' = 0 a new one appears for small mag¬ 
netic helds, where the two dodecahedra are still con¬ 
nected approximately in an antiparallel fashion via the 
J' bonds. This discontinuity increases its strength mono- 
tonically with J', reaching AM ^ 0.8 close to the dimer 
limit (Figs. [3] and |4]). The magnetization curve for 
ij- = J (a = |) is shown in detail in Fig. [S] For small 
fields the spins associated with the non-frustrated bonds 
on the average do not respond as strongly as the rest, 
maintaining a smaller total projection along the z axis 
and having interdodecahedral correlations only weakly 
deviating from -1 (lower right inset of Fig. [S]). The rest of 
the nearest-neighbor correlations deviate more strongly 
from their zero field value, and this deviation increases 
with h. For = 0.095350 the magnetization discon¬ 
tinuity originating in the interdodecahedral coupling ap¬ 
pears. The lowest energy configuration right after the 
discontinuity is similar to the lower-held ground state 
conhguration of an isolated dodecahedron [i^. Even 
though the total magnetization of the two dodecahe¬ 
dra increases right after the low-held discontinuity, the 
net magnetization of the spins connected via J' bonds 
changes direction and points away from the held. These 
spins now share a common polar angle which starts out 
bigger than ^, and then monotonically decreases with the 
held. The hve J' correlations are equal and become more 
antiferromagnetic with increasing held (lower right inset 
of Fig. O. For a specihc value of the held the common 


polar angle becomes equal to and then the J' bonds 
connect antiparallel spins which have zero net magneti¬ 
zation. If the held is further increased the polar angle of 
the J' spins becomes less than ^ and they start to devi¬ 
ate from being antiparallel, while their net magnetization 
is now non-zero and points towards the held. 

Apart from the low-held discontinuity, the net mag¬ 
netization of the J' spins decreases also at the second 
and last ones (Fig. [5] and its upper left inset). The 
three higher-held discontinuities in Fig. [5] are directly re¬ 
lated to the ones of an isolated dodecahedron [la, only 
renormalized by the interdodecahedral interaction. One 
of them, the middle one, is the strongest among all dis¬ 
continuities, and for small J' it is associated with a mag¬ 
netization change AM ~ 1.5 (Figs. OandH]). 

When ^ = 0.84139 (a = 0.45693) the third discon¬ 
tinuity splits up in two (Fig. [5]), which have smaller 
magnitude (Figs. [3]and|4]). The net magnetization of 
the J' spins is shown in Fig. [BKa) for ^ = | (a = |). 
The two discontinuities result in a stepwise increase of 
the magnetization of the J' spins, as seen in the right 
part of the hgure. The rest of the spins also increase 
their net magnetization stepwise. The lower of these two 
discontinuities disappears for ^ = 2.0800 (a = 0.67532) 
(Fig. El). For ^ = 2.78174 (a = 0.73557) and right 
above the second discontinuity two new discontinuities 
emerge (Fig. Eland inset), bringing the total number to 
a maximum equal to six. The net magnetization of the 
J' spins is plotted in Fig. EKb) for the case ^ = 2.82 
(a = 0.73822). Here the net magnetization of the right- 
and left-dodecahedron J' spins is not the same after the 
first new discontinuity, as seen in the right part of the fig¬ 
ure, and the net magnetization of the J' spins decreases 
when this difference appears. The top two discontinuities 
merge for ^ = 3.2552 (a = 0.76499) (Fig. E])) while for 
^ = 3.99890 (a = 0.79996) the second and the third dis¬ 
continuity merge. Finally for ^ ^ 17.1 (a ^ 0.945) the 
second discontinuity changes from a magnetization to a 
susceptibility discontinuity. 

The lowest energy configuration spin directions change 
discontinuously like the magnetization as the gaps are 
encountered with increasing field. Below the first mag¬ 
netization discontinuity the spin configuration is highly 
asymmetric, with each spin having its own polar angle. 
Low symmetry is a general feature of the configurations 
of Fig. El as the spins at best have their own polar angle 
value within an individual dodecahedron, with each of 
these polar angle values shared only by another spin in 
the other dodecahedron. The most symmetric of these 
configurations is when the spins mounted exactly at the 
same location in the two dodecahedra share the polar 
angle, and their azimuthal angles differ by tt. The most 
notable exception to these cases is the lowest energy con¬ 
figuration right after the low-field magnetization jump, 
which is indicated with (blue) up triangles in Fig. El 
which is also the last configuration just before satura¬ 
tion, occuring for fields higher than the ones depicted 
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with (green) diamonds. This configuration is shown in 
Fig. 0 Similarly to the low-field ground state config¬ 
uration of an isolated dodecahedron there are four 
distinct polar angles for the spins, with each one corre¬ 
sponding to a different circle type. Lines of the same type 
represent equal nearest-neighbor correlations. All the az¬ 
imuthal angles are integer multiples of ^, while successive 
azimuthal angles within the same theta group differ by 
Along the central line defined by spins 1, 6, 11, and 
18, nearest-neighbors differ by tt in the azimuthal plane. 
Spins symmetrically placed with respect to this line have 
azimuthal angles adding up to 2tt. The polar angles are 
the same for spins placed exactly at the same locations 
in the two dodecahedra, while their azimuthal angles dif¬ 
fer by TT. In the lowest energy configuration before the 
(orange) x’s in Fig. [5]the polar angles are different in the 
two dodecahedra, and there are 12 distinct polar angles 
for each one of them. 


For s = i an isolated dodecahedron has a discontinu¬ 
ity in the ground state magnetization where its total z 
spin sector = 5 with five spin flip s from saturation is 
never the ground state in a field [2. This results from 
the energy difference of the = 5 and 6 lowest en¬ 
ergy states being smaller than the one of the = 4 
and 5 lowest energy states. The discontinuity carries over 
to the case of the two linked dodecahedra: their lowest 
energy wavefunction for a specific when J' = 0 is 
the product of the individual dodecahedra ground state 
wavefunctions for specific that minimize the en¬ 

ergy and have spins adding up to , with the total en¬ 
ergy equaling the sum of the corresponding individual 
dodecahedra energies. As a result the isolated dodecahe¬ 
dron discontinuity generates two discontinuities for the 
coupled dodecahedra, where the lowest energy levels with 
= 9 and = 11 are never the ground states in a field. 
When J' becomes non-zero, the discontinuities will sur¬ 
vive at least for weak values of it. The ground state mag¬ 
netization curve of Hamiltonian o is calculated for the 
whole range when J' is weak with perturbation the¬ 
ory. The magnetization is also calculated with Lanczos 
diagonalization for arbitrary J', however due to compu¬ 
tational requirements it can not be determined for lower 
in this case. 


A. Perturbation Theory 

For small J' the lowest energies are calculated for ev¬ 
ery sector within first order perturbation theory. The 
unperturbed wavefunction (J' = 0) is the product of the 
lowest energy wavefunctions of the two dodecahedra ac¬ 
cording to the sector they belong to (see Table V 
of Ref. (3). When is away from the single dodec¬ 
ahedron discontinuity for both dodecahedra, the zeroth 


order wavefunction for even is the product of the low¬ 
est energy state with ^ for each dodecahedron: 

= |$^(^))|‘I>^o(y)) (2) 

The index i, j = 1,... ,d counts the degeneracy d of the 
single dodecahedron wavefunction, therefore the unper¬ 
turbed wavefunction of the two dodecahedra is in princi¬ 
ple also degenerate. For odd the combining single 
dodecahedron lowest energy states have 
and ^ ■ Here apart from the degeneracy originating 

in the degeneracy of the single dodecahedron lowest en¬ 
ergy states, an extra factor of two comes about from the 
two distinct values to be accomodated on the two 
dodecahedra. The only exception in the general pattern 
is = 10, where the participating single dodecahedron 
lowest energy states do not have the same = 5 as in 
the other even cases, but due to the single dodecahedron 
discontinuity they have = 4 and 6. Thus in general 
degenerate perturbation theory is required, unless is 
even and different from 10 and the lowest single dodec¬ 
ahedron energy level for ^ is singly degenerate. The 
perturbative term is scaled with J' in Hamiltonian o, 
and first order degenerate perturbation theory produces 
the following matrix, exemplarily for the even case of 
Eq. ©: 

Hi{i * d + j,k * d + 1) = 

5 

m—1 

(3) 

The perturbative term includes combinations of raising 
and lowering operators, as well as diagonal terms. For 
even 10 only diagonal terms can generate non-zero 
energy contributions in first order perturbation theory, 
irrespective of the degeneracy of the single dodecahedron 
state |4 ’q(^)). For odd and = 10 combinations of 
raising and lowering operators also contribute. 

The first order perturbation theory correction for the 
energy Ei is listed for the different sectors in Table 
m According to what has already been mentioned in this 
Subsection about the sectors that combine to form 
the unperturbed wavefunction, away from the disconti¬ 
nuity the energy difference between levels S^, — 1 and 

—1, —2 is the same for J' — 0 when is even. Con¬ 

sequently what determines the relative value of successive 
energy differences between adjacent sectors for weak 
J' are the perturbative energy corrections listed in Table 
m If the relative energy difference between three succes¬ 
sive sectors in Table U when starting from an even 
increases for decreasing S'^, then a new magnetiza¬ 
tion discontinuity appears. This is the case for = 14, 
and the magnetization in a field switches from = 14 
directly to = 12 at least for small J'. This disconti¬ 
nuity does not relate to the one of the isolated dodecahe¬ 
dron but originates in the interdodecahedral coupling. It 
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is then concluded that at least for weak J' and between 
= 8 and 14 the magnetization changes with an exter¬ 
nal field in steps of AS^ = 2. This shows that for two 
linked dodecahedra the magnetization can be changed 
in a controlled way in steps of either AS^ = 1 or 2 by 
adjusting the range of an external magnetic field. 

B. Lanczos Diagonalization 

The magnetization response of Hamiltonian o can be 
calculated for J' of arbitrary strength with Lanczos di¬ 
agonalization, taking into account the spatial sym¬ 
metry of the Hamiltonian (Tbl - fTsL 1^ . In this way 
the Hamiltonian is block-diagonalized according to the 
irreducible representations of its symmetry group. This 
results in eigenstates well-defined according to symme¬ 
try, as well as a Hamiltonian divided in smaller subblocks 
that are easier to diagonalize. In contrast with the per¬ 
turbation theory calculation of Sec. IIV Al here there is no 
restriction on the strength of J', however the calculation 
is limited to higher values of due to computational 
requirements. 

Fig. [5]shows the ground state magnetization curve for 
four different J' values. When J' is small three succes¬ 
sive discontinuities are expected according to Sec. IIV A1 
with the subsectors S'^ = 9, II and 13 never including the 
lowest energy state in a magnetic field. Two discontinu¬ 
ities are highlighted with (red) arrows in Fig. EKa) where 
J' = while the one that corresponds to the lowest 
is not mapped out due to the computational require¬ 
ments to hnd the lowest energy state for = 8. When 
J' = ^ (Fig. EKb)) the interdodecahedral coupling is not 
weak any more, and the highest discontinuity has dis¬ 
appeared. For even higher J' — ^ (Fig. EKc)) the dodec¬ 
ahedra feel each other’s influence more strongly, and as a 
result the intermediate discontinuity disappears as well. 
Fig. EKd) shows the magnetization for J' = J where all 
couplings are equal and even magnetization plateaus are 
absent for the range of the weaker J' discontinuities. 

Fig. E] shows the distinct ground state expectation 
values of the nearest-neighbor correlation functions < 
Si ■ Sj > for J' = -^ and There are in principle six 
unique such correlations, with one for each of the rings 
that respectively contain spins I to 5, 6 to 15, and 16 to 20 
(Fig. [7]), and two more for nearest-neighbor correlations 
between spins that belong to different rings. Tha last 
unique correlation is between spins belonging to differ¬ 
ent dodecahedra. For J' = -^ the single-dodecahedron 
character is preserved, at least not for high 5^, where 
only correlations within the two dodecahedra are anti¬ 
ferromagnetic. The intradodecahedral correlation (rep¬ 
resented by < Si • 5*21 > in Fig. El) is ferromagnetic 
for all the presented. The situation is different when 
J' = ^. Now the intradodecahedral correlation acquires 


an AFM character, which is the strongest along with its 
neighboring correlation < Si ■ S 2 >■ 

Fig. (TU] shows the distinct ground state expectation 
values of the projections of the individual spins along 
the z axis < sf > for J' = -^ and There are in 
principle four unique such projections, with spins I to 5 
having a common one (Fig. [7]), and spins 16 to 20 an¬ 
other one. Also every second of spins 6 to 15 shares the 
same value of < sf >. Spins of the central and the outer 
pentagon have lower values for J' = which agrees 
with their stronger intradodecahedral AFM correlations 
of Fig. EKa)- For J' = ^ the central pentagon has even 
lower < sf >, which now corresponds to the strongest 
AFM correlations being between spins in this pentagon, 
and between spins in this pentagon and their counter¬ 
parts in the corresponding pentagon of the other dodec¬ 
ahedron connected via the J' bonds, as shown in Fig. 

Elb). 


V. CONCLUSIONS 

The ground state magnetic response of two coupled do¬ 
decahedra was investigated within the framework of the 
AHM. The classical magnetization discontinuities of an 
isolated dodecahedron were found to be renormalized by 
the interdodecahedral coupling, while new ones emerge. 
For a specific range of the coupling the total number of 
discontinuities goes up to six. At the full quantum limit 
s = b the isolated dodecahedron magnetization disconti¬ 
nuity gives rise to two neighboring discontinuities, with 
a third one appearing adjacent to these two. The two 
dodecahedra system has a magnetic response which for 
a significant range of the field is associated with magne¬ 
tization steps with AS^ = 2, which is twice as strong as 
the usual magnetization difference between adjacent 
sectors for a quantum spin system. This shows that the 
magnetization change can be controlled by adjusting the 
range of an external magnetic field. 

The frustrated nature of the dodecahedron results in 
unexpected ground state magnetization discontinuities in 
a held when the AHM is considered on it. Usually such 
discontinuities are associated with magnetic anisotropy, 
but in this case they are allowed by the special connec¬ 
tivity of the dodecahedron. The formation of a two do¬ 
decahedra molecule with the introduction of unfrustrated 
coupling between one of their faces enriches the ground 
state magnetic response. Along these lines, it is of in¬ 
terest to extend this investigation on more than two do¬ 
decahedra linked together to form a chain-type or even 
more complicated structures, and calculate the magnetic 
response while the individual cluster frustration and the 
coupling between clusters compete in the presence of an 
external magnetic field. 
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TABLE I: First order perturbation theory energy correction 
El for the different sectors for the two coupled dodeca- 
hedra. The numbers were generated with double precision 
but are presented with five significants digits for the sake of 
brevity where applicable. 


S" 

El 

S" 

El 

S" 

El 

0 

0 

7 

0.0045106 

14 

0.51917 

1 

-0.15988 

8 

1 

5 

15 

0.47352 

2 

0.0052328 

9 

0.11821 

16 

0.62767 

3 

-0.064727 

10 

3 

in 

17 

0.49119 

4 

-0.083283 

11 

0.32070 

18 

0.88820 

5 

-0.16495 

12 

9 

20 

19 

0.88820 

6 

-0.0039974 

13 

0.48578 

20 

5 

4 


FIG. 1: Projection of a single dodecahedron on a plane. The 
solid lines are antiferromagnetic bonds with strength J and 
the circles spins with magnitude s. The cluster considered in 
this paper consists of two such dodecahedra connected with 
J' interactions between the spins of one of their faces, high¬ 
lighted with white color and including spins 1 to 5. The spins 
of the second dodecahedron have the same linking pattern 
with the first as shown in the figure and their indices are 
given as 20 -I- i, where i runs from 1 to 20. The J' bonds are 
between spins i and 20 -I- i, with i running from 1 to 5. 




FIG. 2: (Color online) Location of the classical ground state 
discontinuities with respect to a and . Each discontinuity 
is distinguished by a different symbol (and color). The (black) 
circles, (red) squares and (green) diamonds correspond to the 
magnetization discontinuities that originate in the isolated 
dodecahedron (J' = 0). The (maroon) circles on the far right 
correspond to the sole susceptibility discontinuity. The inset 
shows in detail the three magnetization discontinuities that 
are close to each other. 



FIG. 3: (Color online) Inaccessible lowest energy state mag¬ 
netizations per spin ^ as functions of a. The inaccessible 
magnetizations are distinguished by their color, in accordance 
with Fig. m 
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FIG. 4: (Color online) Magnetization change per spin 
for the classical ground state magnetization discontinuities 
as a function of a. Each discontinuity is distinguished by a 
different symbol (and color), in accordance with Figs. [5] and 

11 



FIG. 5: (Color online) Total classical ground state magneti¬ 
zation projection along the « axis per spin ^ as a function 
of for ^ = y (a = |). The (black) continuous line cor¬ 
responds to the spins in each dodecahedron connected via J' 
bonds, the (red) long-dashed line to the rest of the spins, and 
the (green) long-dashed-dot line to all the spins. The lower 
right inset shows the average interdodecahedral correlation 
function | < s; • Sj > |j/ for small fields. The upper left inset 
shows in detail two of the discontinuities. 
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FIG. 6: Total classical ground state magnetization projection 
along the z axis per spin ^ as a function of for the spins 

in each dodecahedron connected via J' bonds for (a) -^ = | 
(a = |), and (b) ^ = 2.82 {a = 0.73822). The long-dashed 
line gives the average when the two values differ. 



FIG. 7: (Golor online) Projection of a single dodecahedron 
on a plane (see Fig. [T] for description). Lines of the same 
type (and color) represent equal nearest-neighbor correlations 
Si ■ Sj in the lowest-energy configuration right after the low- 
field magnetization discontinuity and just before saturation 
for classical spins, and equal nearest-neighbor correlations 
< Si ■ Sj > in the ground state for s = | for different 5^. 
Circles of the same pattern (and color) represent equal polar 
angles in the lowest-energy configuration right after the low- 
field magnetization discontinuity and just before saturation 
for classical spins, and equal projections of local spins along 
the z axis < sf > in the ground state for s = ^ for different 
5L 
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FIG. 8: (Color online) in the ground state as a function of 
for s = I and J' equal to (a) (b) (c) and (d) J. 


The (red) arrows point at the magnetization discontinuities, 
where AS^ = 2. 





FIG. 9: (Color online) Distinct ground state expectation val¬ 
ues of the nearest-neighbor correlation functions < Si ■ Sj > as 
a function of for s = | and J' equal to (a) and (b) 

The spin indices of the legend (and the color coding) refer to 

Fig.m 
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FIG. 10: (Color online) Distinct ground state expectation 
values of the projections of local spins along the 2 axis < s| > 
as a function of for s = ^ and and J' equal to (a) and 
(b) The spin indices of the legend (and the color coding) 
refer to Fig. [7] 




